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In this paper we investigate the general features of " Oscillatory Inflation" . In adiabatic approxi- 
mation , we derive a general formula for the number of e-foldings N which reduces to the standard 
expression in case of the slow role approximation and leads to the Damour-Mukhanov type expres- 
sion for the slowly varying adiabatic index. We apply our result to the logarithmic potential and 
arrive at a simple and elegant formula for the number of e-foldings. We evolve the field equations 
numerically and observe a remarkable agreement with the analytical result. 



I. INTRODUCTION 

o : 

The inflationary universe scenario has become an integral part of the standard model of universe. It does resolve 
(N ■ the outstanding problems of the standard model like the flatness problem, the horizon problem, the homogeneity and 
isotropy problem etc. It also provides an important clue for the origin of structure formation in the universe [1,2, 3] . 
The underlying idea of inflation is that there was an epoch when the universe was slowly rolling down the flat wings 
of the scalar potential such that vacuum energy was dominant leading to the exponential growth of the scale factor. 
At the end of slow roll, the universe falls into the core of the potential and fast oscillates near the minimum leading 
to the re-heating ultimately. Damour and Mukhanov have discussed the possibility of inflation during the period of 
oscillations [ 4 ] . Further discussion on the same theme can be found in reference [ 5 ]. In fact, such a possibility was 
first pointed out by Turner[ 6 ]. Inflation with oscillation can be realized by a non-convex potential V(tf>) ( < 0) 

■ with small convex core near the minimum of the potential and with large flat wings such that universe spends most 
[--««. ' of the time away from the core and inflates. In this paper we study general aspects of inflation during oscillations . 

. We derive general expression for the number of e-folds which reduce to the standard formula in the slow role regime 
and leads to the Damour-Mukhanov type expression for slowly varying adiabatic index. We apply our result to the 
logarithmic potential and show that the model exhibits simple analytical solution which has remarkable agreement 

• with the numerical simulation. 

"o 

II. EVOLUTION EQUATIONS AND CRITERIA OF INFLATION WITH OSCILLATIONS: 

U '. 

The field evolution equations in the Freedman cosmology have the form 

H 2 = -%p (2) 



a 3M 2 
where 



= -3H(p + p) =-3H0 2 (3) 

a 47T , 

(p + 3p) (4) 



9 = \fi + V(4>) (5) 



P = \tf - v(4>) (6) 

Equations (3) and (4) are not independent of (1) and (2) but will be useful in view of forthcoming discussion. Equation 

52 

(3) reveals that the energy density in scalar field 4> decreases because of the red shifting away of the kinetic part . 
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In what follows , we shall assume that the potential is even and has minimum at (f> — . When the field initially being 
displaced from the minimum of the potential, rolls below its slow roll value, the coherence oscillation regime u > > H 
commences. The evolution equation can then be approximately solved by separating the two times scales namely the 
fast oscillation time scale and the longer expansion time scale. On the first time scale, the Hubble expansion can be 
neglected, and one obtain ^ as a function of time; 

t-t = ± [ 1 =#, (7) 

J y/2(V m - VW) 

where p = V m = V((f> m ) ; V m being the maximum current value of the potential energy and (j) m being the field 
amplitude . On the longer time scale p and 4> m slowly decrease because Hubble damping term in equation ( 1 ) .The 
average adiabatic index 7 is defined as , 

In I m\ I J.2 \ 

(8) 

where< . > denotes the time average over one oscillation . Equation (4) then tells that expansion during oscillations 
would continue (a > 0) if 7 < |. The adiabatic evolution of a{t) and p is given by , 

a(t) cx t& (9) 




p = V{4> m ) cx tr 2 (10) 
As <j> — ^p,the condition 7 < | can equivalently be written [ 6 ] , 

,_/£\_<^>_ J(1 _<v> 

\ P I V m V m 

f* m {i-v{<f>)/v m )-*d4> 3 

It is suggestive to write (11) in an equivalent form, 

{U(4>)) = (V - > (12) 

The expression (12) has a simple geometric meaning; the average over one oscillation of the intercept U(</>) at <f) should 
be positive to ensure accelerated expansion. In models of inflation with oscillation the slow roll value of the field 4> s is 
" far away" from convex core given by <j) c such that(C/) is positive before (f> rolls below certain critical value of the field 
4> cr Fig. 1 . As a result the universe undergoes accelerated expansion (a > 0). This is because the average of intercept 
U((f>) is negative for < < (j) c and the positive contribution to < U(<p) > starts coming only as <f> gets bigger than 
<p c . Obviously one has to take <j> further to some <f> C r in order to cancel the negative contribution to < U (4>) > coming 
from the convex core . In certain sense the <f> cr defines the effective size of convex core . Clearly , as <f> rolls below 
4> cr , the average of U{4>) over one oscillation turns negative and consequently inflation ceases loosing to oscillations 
without accelerated expansion . Hence for inflation to occur during oscillations , the field amplitude 4> m should very 
between <fi cr and <j> s . We shall call this regime " oscillatory inflation regime " . 
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FIG. 1. Qualitative picture of an effective potential showing the intercept U(<j)) of the tangent to the curve V(<j)) at various 
values of the field 0.For <f> < </) c ,the value of the intercept is negative . U(4>)is positive for <f> > <f> c . The critical field 4> cr is fixed 
such that the positive contribution coming from the interval <f> c < 4> < 4> cr to the average intercept (U(4>)) exactly cancel the 
negative contribution of the convex core . {U(4>)) is positive as <f> varies between <f> cr and 4> s and the universe expands in this 
regime called the " oscillatory inflation regime " . 



III. MODEL: 

We shall study here the logarithmic potential which may implement inflation with oscillations, 

V(4) = ^oln[l + (^ c ) 2 )] 
This is the q tending to zero limit of the potential suggested by Damour and Mukhanov, 



q/2 



(13) 



(14) 



The potential has a convex core near the bottom of potential given by 4> c whereas is negative away from the 
core. For >> 4> c ,V(4>) = 2Vq ln(4>/<j> c ) and consequently 7 ~ l/ln(0 m /0 c ) = l/ln(/3) 



A. critical value of field 



In order to compute the size of inflation , it would be necessary to determine (U) : 

(U) = [\v(x,(3)-xV x (x,(3))dx 
Jo 



(15) 



where (5 = (f> m /4>c , x = <p/4> m and V(x, [}) — Vo ln(l + /3 2 x 2 ) in case of the potential given by ( 13 ). Since only the 
sigh of (U) is important we have omitted several positive constants in the expression( 15 ) . The average intercept 
(U) is plotted in Fig. 2 in case of the logarithmic potential under consideration .The value of (3 where (U) changes 
sign precisely determines the critical value of the field or equivalently /3 cr .For the logarithmic potential /3 cr turns out 
to be nearly equal to 3.2. 
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FIG. 2. Plot of the average intercept (U(<f>)) versus (3 = 4> m /4>c for the potential given by Eq.( 13 ).{U ((/>)) changes sign as 
m becomes nearly equal to 3.20 c indicating the beginning of oscillatory inflation regime. 

B. slow roll field: 

The slow roll parameter e is defined as : 



_ AP 
6 ~ 16tt 
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(16) 



V 

The slow roll ends when e ~ 1. Using ( 16 ) yields the slow roll value of the field cf) s , 

Vl6^x (3 s \nf3 s ~M/q> c (17) 

where (3 S = ^f- and (f> s depends upon <p c .It should be emphasized that (3 S or equivalently <p s would always depend 
upon 4> c in case of non-power law potentials. 

C. Number of e-folds: 

The number of e-folds is defined by , 

7V = ln^, (18) 

Since the scale factor "a" depends upon (j> m or equivalently upon (3, ( 18 ) should be used carefully . We divide the 
interval {(3 cr ,f3 s ] into small slices , A(3j — f3j+\ — /3j-with j — l,M such that (3\ = [3 cr and (3m — /3 S -The number of 
e-foldings obtained for the interval A/3, is given by , 

where a J+ i < aj for universe rolling down from (3 s to(3cr using ( 9 ) and ( 10 ) leads to the following expression , 

N^p 3 = ^\n(V{P 1+1 )/V{p 3 )) 
Expanding V((3j + i) in Teller series in the neighborhood of (3j and using ln(l + x) ~ x for x « 1 we get , 
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Integrating ( 19 ) from fi cr to /3 S one obtains the total number of c-foldings , 



= / 3 N(P)d0 



N 

I P. 
where 

V'(0)/V(/3) 



N(0) 



37(/3) 



N(P) can be interpreted as the density of e- foldings with regard to (3. The correct definition , however , should take 
into account the change in the comoving Hubble length [ 7 ] , 



Similar arguments as above can be used to obtain the expression forTV , 

re 

10, 

where 



N= N((3)dp (21) 

JPr.r- 



Expression ( 22 ) formally written for oscillatory regime presents a general result for the number e-foldings . In fact 
,in slow roll approximation , 

7 ~ + V + ~ V ~ \24ir) \ V 

Substituting this expression in ( 22 ) leads to the standard formula for the number of e-foldings [ 2 ]. Needless to say 
that N nearly coincides with N in case of the slow role regime .In case 7 is a slowly varying function of /3, N assumes 
the following form , 



N =(—--) In 



(23) 



For a power law potential ( 23 ) readily reduces to the expression obtained by Liddle and Majumdar [ 6 ]. 

IV. NUMERICAL SIMULATION AND ANALYTICAL RESULTS 

We will show that a simple and elegant analytical formula for the number of e-folding can be obtained in case of 
the logarithmic potential ( 13 ). As mentioned earlier that 7 ~ l/ln(/3) in the present case [ 4,5 ]. Substituting the 
approximate value of 7 in ( 22 ) leads to a simple expression for N(/3), 

*w = h-2itp ^ 

Using ( 21 ) we obtain the formula for the total number of e-foldings , 

3 /3 cr 2 \ln[3 cr ) V ' 
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It should be mentioned that we have used the approximate value of 7 to obtain the above expression for TV. Since 
the approximation holds good away from the core and little contribution to the number of e-foldings comes from the 
regions near the core, we would expect that the expression ( 25 ) gives the correct value of TV. We have evolved the field 
equations numerically and studied the evolution of the number of c-foldings. Fig2 displays the numerical evolution 
of TV for two different values of c .Thc dependence of fi s on <f> c can in general be red off from ( 17 ).Thc heights of 
maxima of TV are described by ( 25 ).We see the remarkable agreement of our analytical result with the numerical 
simulation. As emphasized earlier [ 4, 5 ] the number of e-foldings TV is small. Even if <j> c is pushed to electro-weak 
scale TV does not exceed 10. One might think that brane world inflation would be helpful here as the brane corrections 
enhances the prospects of inflation [ 8, 9, 10]. Unfortunately , this is not so as the brane induced term increases the 
Hubble expansion rate which is not relevant during oscillations. 
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FIG. 3. Plot of TV versus r\ = tM based upon numerical simulation for two choices of 4> c .The evolution takes place from the 
end of slow roll and counts the number of e-foldings during the oscillations. The heights of maxima give the actual number of 
e-foldings described by ( 25 ). 

We have described the general features of oscillatory inflation. We have been able to obtain a general expression for 
the number of e-foldings and demonstrated its consistency. Applying the same to the logarithmic potential we obtained 
an elegant and simple analytical formula for the number of e-foldings. We evolve the field equations numerically and 
observe a remarkable agreement of the analytical result with the numerical simulation. 
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